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ABSTRACT 

The HERschel Inventory of The Agents of Galaxy Evolution (HERITAGE) of the Magellanic Clouds will use dust emission to investigate the 
<^ ° life cycle of matter in both the Large and Small Magellanic Clouds (LMC and SMC). Using the Herschel Space Observatory's PACS and SPIRE 

O photometry cameras, we imaged a 2° x 8° strip through the LMC, at a position angle of ~22.5° as part of the science demonstration phase of the 
Herschel mission. We present the data in all 5 Herschel bands: PACS 100 and 160 jdm and SPIRE 250, 350 and 500 jdm. We present two dust 
^ ^ models that both adequately fit the spectral energy distribution for the entire strip and both reveal that the SPIRE 500 jim emission is in excess of 
^ the models by ~6 to 17%. The SPIRE emission follows the distribution of the dust mass, which is derived from the model. The PAH-to-dust mass 
I if pah) image of the strip reveals a possible enhancement in the LMC bar in agreement with previous work. We compare the gas mass distribution 
O derived from the HI 21 cm and CO J=l-0 line emission maps to the dust mass map from the models and derive gas-to-dust mass ratios (GDRs). 
^ The dust model, which uses the standard graphite and silicate optical properties for Galactic dust, has a very low GDR = 65^}gmaking it an 
C/) , unrealistic dust model for the LMC. Our second dust model, which uses amorphous carbon instead of graphite, has a flatter emissivity index in the 
submillimeter and results in a GDR = 287^^2 th^t is more consistent with a GDR inferred from extinction. 

. Key words. Galaxies: Magellanic Clouds, ISM: dust, Submillimeter: galaxies, ISM 

> : 

■ 1 ■ Introduction emission from the ISM and stars, which are the agents of galaxy 

evolution. Herschel SPIRE and PACS images provide key in- 
^ . The Large Magellanic Cloud (LMC) and the Small Magellanic eights into the life cycle of galaxies because the far-infrared and 
. ■ Cloud (SMC) are the best astrophysical laboratories to study the submillimeter emission from dust grains is an effective tracer of 
^ lifecycle of the interstellar medium (ISM), because their prox- ^^e ISM dust, the most deeply embedded young stellar objects, 
^ imity (50 kpc, e.g. Schaefer 2008; 61 kpc, Szewczyk et al. 2009) ^nd the dust ejected by evolved massive stars and supernovae. 
^ permits detailed studies of resolved ISM clouds and their re- j^^^^^g ^^^^^^^ demonstration phase (SDP), we imaged 
„ lation to stellar populations on global scales, m an unambigu- ^ ^^^^ ^^^^^ ^^^^^^ LMC that covers the entire range of inter- 
^ > . ous manner, and as a controlled function of environment. Their ^^^^^^ ^^-^^^^ ^^p^^^ ^^^^^ ^.^j^ HERITAGE: giant molec- 
^ sub-solar metallicities (Zlmc ^ 0.5 x Z©, Zsmc ^ 0.2 x Z©; ^^^^ p^j ^^^^^^ (Roman-Duval et al. 2010; Kim et 
r\ Dufour et al. 1982, Bernard et al. 2006) permit investigations 2OIO), HII regions (Hony et al. 2010), young stellar objects 
^ . on how processes governing galaxy evolution depend on metal- ^^^^^j^ ^1 2010), supernova remnants (Otsuka et al. 2010), 
hcity. The Herschel Observatory (Pilbratt et al. 2010) open- ^^^^^^^ ^^^^^ ^B^^^^ ^1 2010). In this paper, we present 
time key program, entitled HERschel Inventory of The Agents ^^e observing strategy, processing and resulting data set (§ 2), 
of Galaxy Evolution (HERITAGE) in the Magellanic Clouds, ^^^^^^^^ ^^^^^^ distribution (SED) of the strip (§ 3) and the 
will perform a uniform survey of the LMC (8° x 8.5°), SMC ^^^^^ ^^^1 distribution of the dust and gas (§ 4). 
(5° X 5°), and the Magellanic Bridge (4° X 3°) with the Spectral ^ ^ & / 
and Photometric Imaging Receiver (SPIRE) at 250, 350, and 
500 yum (Griffin et al. 2010) and with the Photodetector Array 
Camera and Spectrometer (PACS) at 100 and 160 yum (Poglitsch 
et al. 2010). The HERITAGE science goals are to study the life We observed a 2° x 8° strip through the LMC, at a position 
cycle of matter in the Magellanic Clouds by probing the dust angle of -22.5° using the Herschel Observatory instruments 
SPIRE and PACS in parallel observing mode. Observations be- 



2. Observations and data reduction 



* Herschel is an ESA space observatory with science instruments gan on November 22, 2009 at -23:00 h UT and lasted -18 hours, 
provided by European-led Principal Investigator consortia and with im- Two 9 hour astronomical observation requests (AORs) were con- 
portant participation from NASA. structed back-to-back in time to cover the region. The observed 
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Fig. 1. The HERITAGE LMC SDP strip in the five Herschel bands with color bars on top. The strip has been reprojected to orient 
the figures up down, and represent the SPIRE coverage. The RA and Dec gridhnes for all images here and in Fig. 3 are shown on 
the SPIRE 500 fim image. The PACS images are shifted because their coverage is shifted from the SPIRE data in parallel mode. 
From left to right the color scale ranges and tick steps for the arcsinh stretch (Lupton et al. 1999) in units of MJy/sr: to 31, steps 
of 5 (SPIRE 500 yum), to 71, steps of 10 (SPIRE 350 yum), to 161, steps of 20 (SPIRE 250 yum), to 250, steps of 25 (PACS 160 
yum) and -20 to 150, steps of 15 (PACS 100 yum). Last on the right, a 3 color image with SPIRE 350 yum in red, MIPS 70 yum in green 
and MIPS 24 yum in blue (Meixner et al. 2006). 



wavelength bands include PACS 100 and 160 yum and SPIRE 
250, 350 and 500 yum (Fig.ffl). 

2.1. PACS data reduction 

The PACS data were reduced starting from the level product 
using the HIP^ version 2.0 data reduction software (Ott 2010). 
In particular, bad and saturated pixels were masked, and a flat 
field correction and the photometric calibration were applied. 
The baseline signal level of each bolometer was estimated at the 
beginning and end points of each scan leg which lie outside of 
the LMC. For each pixel and each scan leg, a linear fit linking 
these two points was subtracted from the corresponding time- 
line. The multi-resolution median transforms (MMT) deglitch- 
ing algorithm was then applied to correct sudden jumps in the 
timeline, or glitches, due to cosmic ray hits on the detector. 

Initial two-dimensional maps were created for the 100 and 
160 yum bands. However, severe striping appeared along the scan 
direction of the PACS maps. To mitigate the stripes, we sub- 
tracted an image of the striping, produced by unambiguously 

^ RIFE is a joint development by the Herschel Science Ground 
Segment Consortium, consisting of ESA, the NASA Herschel Science 
Center, and the HIFI, PACS and SPIRE consortia. 



identifying the power spectrum associated with the stripes, from 
the PACS images. The peak-to-peak variation was reduced from 
0.021 to 0.01 1 in the PACS 100 and 160 yum images (Fig.[T]). We 
do not use the PACS data for analysis in this paper, but antici- 
pate that with cross-scans the PACS data will be viable for an 
extended map of the diff'use ISM. Nevertheless, the de-striping 
was eff'ective for small regions (Otsuka et al. 2010) and revealing 
point sources (Sewilo et al. 2010; Boyer et al. 2010). 

The PACS 160 yum brightness values are systematically 
15% lower than the corresponding MIPS 160 yum values from 
[Meixner et al. (2006)| using the adopted HIPE 1.2 calibration 
tree which is consistent with the 20% absolute flux calibration 
error for PACS (Poglitsch et al. 2010). The accuracy of the PACS 
astrometry was estimated by comparing the positions of point 
sources detected in Spitzer and PACS images of the LMC. We 
find off'sets up to 6.5'', with no preferred direction, which have 
been traced to issues with the star trackers during very large scan 
maps. 

2.2. SPIRE data reduction 

The SPIRE data were processed with HIPE version 2 (Ott 2010) 
and updated calibration products provided by the SPIRE team. 
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Fig. 2. The spectral energy distribution of the LMC HERITAGE 
SDP strip. Open diamonds with error bars mark the fluxes from 
the SPIRE 500, 350, and 250 fim, and the SAGE-LMC (Meixner 
et al. 2006) MIPS 160, 70 and 24 jim and the IRAC 8, 5.8, 4.5 
and 3.6 fim bands in solar luminosities. The grey solid line is the 
result of model 2 using the amorphous carbon dust (see § 3). The 
green filled points are the simulated photometry based on model 
2. The ISM dust emission components are shown in blue, green 
and red according to the legend. The yellow line is the stellar 
emission component. 



Using custom routines, we removed jumps in the timeline caused 
by co-occurring glitches that aff'ect all detectors of a single array 
simultaneously and by jumps in the thermistor voltages. We in- 
clude the data at the end of each raster leg, during which the tele- 
scope turns around, in order to increase the coverage. For each 
scan, for each detector, we derived the median, baseline value 
of the measured flux at each end of the scan and subtract a lin- 
ear baseline fitted to the two end points from the scans. Images 
at 250, 350 and 500 yum were constructed using the same refer- 
ence point near the map center and are shown in Fig.[TJ In order 
to remove any instrumental residuals in the SPIRE images, we 
use regions at the ends of the map, assume they should have 
zero emission and subtract a linear gradient from the images. 
This process was done in a consistent way for all ancillary data 
(MIPS, IRAC images and HI 21 cm emission images) so that 
comparisons between the data sets during the analysis are done 
in a consistent fashion. 

As recommended by the SPIRE ICC we multiplied the maps 
by the flux calibration correction factors 1.02, 1.05, 0.94 for 
250, 350 and 500 yum, respectively. The uncertainty in the fi- 
nal absolute point source flux calibration is ±15% (Swinyard et 
al. 2010). The final SPIRE maps were converted from Jy per 
beam to MJy/sr using the eff'ective beam areas of 9.28x10"^, 
1.74xl0-^ and 3.57x10"^ sr for 250, 350 and 500 yum. The po- 
sitions of point sources in the 250 yum map were compared with 
their positions at 24 yum from MIPS revealing off'sets up to 10'', 
similar to PACS. We corrected this problem with a small spheri- 
cal rotation of the map, although some small residual astrometry 
errors remained, particularly far from the map centers. 

3. Spectral energy distribution of the strip 

The SED of the whole strip (Fig. O is fit with a com- 
prehensive dust radiation model described in detail by 
[Galametz et al. (2009), We adopt the Dale & Helou (2002; Eq. 1) 
prescription, relating the dust mass to the integrated energy den- 
sity to which it is exposed in units of the solar neighborhood 
intensity (2.2 x 10"^Wm-2; Mathis et al. 1983). The standard 



model (model 1) we apply on the strip uses the Galactic ISM dust 
size distribution and composition of graphite and silicate grains 
from Zubko et al. (2004; BARE-GR-S model) who successfully 
fit the Galactic ISM dust emission out to submillimeter wave- 
lengths. The dust emissivity index at the longest wavelengths is 
/3 = 2, which is also t he assumed emissivi ty value for the SAGE- 
LMC ISM study by Bernard et al. (2008), However, as wiU be 
shown in Sect. 4, the grains of this model do not have enough 
emissivity in the submillimeter regime to fit the SPIRE fluxes, 
with a gas-to-dust mass ratio (GDR) consistent with the elemen- 
tal abundances of the LMC. As a consequence, we also created 
another model (model 2) replacing the graphites (of model 1) 
by amorphous carbons from Rouleau & Martin (1991; ACl). 
This modified model has more emissivity at long wavelengths 
and therefore gives a realistic GDR. For both models, the PAH- 
to-dust mass fraction, /pah is reported in units of the Galactic 
PAH fraction of 0.046 (Draine & Li 2007); i.e. /pah = 1 means 
Galactic PAH abundance. 

Both dust models provide an adequate fit to the SED of the 
strip. For model 1, the PAH-to-dust mass fraction is slightly 
lower than the Galactic mass fraction (/pah = 0.85^^ ^^), and its 
mass average stellar light intensity is twice the solar neighbor- 
hood intensity ({U} = 2.0^^^). Although Hii regions can show 
much hotter dust (Hony et al. 2010), it does not dominate the 
emission. We had to decrease the mass fraction of small grains 
(< 10 nm) by a factor of 2, in order to fit the 24 yum flux of the 
difl'use ISM. This has no significant eff'ect on the total dust mass. 
Model 2 gives a lower PAH mass fraction {/pah — ^-^^^J^qq^) and 
higher starlight intensity ({U} = 9.1!|-^). Notably, the 500 yum 
flux is in excess compared to our SED model of the whole strip 
by -17% (model 1) and ~6% (model 2). The origin of this 
excess is still unknown; however, such submillimeter excesses 
have been observed before in other magellanic irregular galaxies 
by e.g. Ga lliano et al. (2005 ), The spatial distribution and poten- 
tial locati on of this 500 fim flux excess is discussed in further 
details by |Gordon et al. (2010)t 



4. Distributions of dust and gas in the strip 

The emission at these Herschel wavelengths (Fig. [T]) is domi- 
nated by dust emission from the ISM and similar to the Spitzer 
MIPS 160 yum emission (see e.g. Bernard et al. 2008). In fact, 
the SPIRE maps have an identical appearance in the overall strip 
which is expected because they trace the Rayleigh- Jeans tail of 
the ISM dust emission (Fig.[2l). The three color image in Fig.[T] 
shows the variation in color temperature of this ISM dust emis- 
sion. The SPIRE 350 yum (red) samples all dust emission, includ- 
ing some of the coldest dust (T ^ 10 K), the MIPS 70 yum (green) 
is near the SED peak, and the MIPS 24 yum (blue) represents the 
hottest dust emission peaking in the H II regions. 

The modelling described in Sect. 3 has also been applied 
to each pixel of the strip (84'' in angular size or 20 pc in lin- 
ear size). The dust mass surface density, ^dust^ average stellar 
light intensity, (t/), and PAH-to-dust mass ratio, fpAH, derived 
from dust model 2 are shown in Fig. [3j The details of these 
models are described by Galliano et al. (in prep.). We also de- 
rive a gas mass surface density (J^^gs ) image using the HI 21 
cm emission from ATCA-hParkes by |Kim et al. (2003)[ which 
traces the atomic neutral gas, and the CO J=l-0 emission from 
the NANTEN survey by |Fukui et al. (2008)[ which traces the 
molecular gas. The HI (T resolution) and CO (2.6' resolution) 
maps were convolved to the IRIS 100 yum resolution (4.3') using 
Gaussian kernels of widths 4.18' and 3.42' respectively. We used 
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Xco = 7 X lO^^cm-^K-^kms-i from Fukui et al. (2008) to de- 
rive the H2 mass from CO. We combine the HI and H2 column 
densities and adjust for the contribution of heHum (xl.36 more) 
to create the total gas mass map in Fig. 3. 

The dust mass distribution is spatially coincident with 
the SPIRE emission (Fig. [3]) and is discussed in detail by 
[Gordon et al. (2010)1 The HI 21 cm emission appears almost 
identical to the dust mass i mage (Fig. [3]) and the similarity is 
quantified in more detail by |Kim et al. (2010)| The CO J=l - 
emission coincides with the peaks in the gas mass (Fig. O 
and t wo molecular clou ds are discussed in more detail by 
[Roman-Duval et al. (2010), The average starlight intensity heat- 
ing the dust, (U), approximately follows the IRAC 3.6 jim image 
(e.g. Meixner et al. 2006) being enhanced in the stellar bar, and 
in the giant HII region complexes, such as N44, which is just 
north of the stellar bar. The spatial distribution of /pah derived 
by the model appears to be enhanced toward the stellar bar rela- 
tive to the rest of the ISM (Fig. 3) in agreement with the findings 
of [Paradis et al. (2009)[ We calculate that this enhancement is 
equivalent to < 2000 M© (model 1) or < 300 M© (model 2) in 
PAH. However, it could be the result of the increased hardness of 
the starlight along the bar, mistaken for an increased abundance. 

Over the area in common between the gas and dust mass 
images, we calculate the total gas mass and divide by the to- 
tal dust mass to derive a GDR. For model 1 , we calculate a 
GDR = 65^} ^. This value is less than the Galactic GDR value 
of 157 from [Zubko et al. (2004)| which would be unexpected 
given the LMC's lower metallicity. In fact, the model 1 GDR 
is a factor of 3 low er than the SAGE-LMC ISM study by 
[Bernard et al. (2008)} who adopted similar Galactic dust proper- 
ties. It is also too low by a factor of three compared to UV extinc- 
tion measurements of |Gordon et al. (2003)[ Thus, the addition of 
the SPIRE data demands too much dust, if we assume standard 
graphite and silicate optical properties which have been success- 
fully used for Galactic dust models. For model 2, we calculate a 
GDR = 287^42' which is much more consistent with prior obser- 
vations. Essentially the required dust mass for model 2 is smaller 
by a factor ~3 because it requires a larger fraction of the dust to 
be illuminated by intenser radiation fields and therefore less dust 
mass to account for the observed SPIRE measurements. The fact 
that the pixel to pixel averaged intensity (U) = 7.6^3 ^ for model 
2 is higher than for model I ({U} = 1.0^^4) makes model 2 more 
consistent with the expected intense radiation conditions within 
the LMC. The amorphous carbon emissivity is flatter in the sub- 
millimeter (fi < 2 ) which is in agreeme nt with the independent 
approach taken by [Gordon et al. (2010)| for this data set and with 
the independent analysis of the TOPHAT and DIRBE measure- 
ments of the LMC by |Aguirre et al. (2003)[ 

Although our results indicate that the standard graphite and 
silicate optical properties for Galactic dust are not appropriate 
for the LMC dust, our suggested amorphous carbon and silicate 
dust model is not necessarily unique or the best model. The anal- 
ysis of dust models begun in this paper will need to be revisited 
when the full HERITAGE data set is available. In particular a 
full exploration of parameter space for dust properties including 
composition and size distribution can be investigated in terms of 
ISM environment and metallicity. 
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Fig. 3. Analysis images of the SDP strip based on model 2 
(amorphous carbon, see § 4) with color bars for each on top, 
described from left to right. Dust surface mass, Y^dust^ in red 
color arcsinh scale ranging from to 0.255, in tick steps of 0.05 
M© pc"^ with HI contours at levels 1, 2 and 4 xlO^^ H cm"^ 
from |Kim et al. (2003)| Gas surface mass, S^^^, in jet color arc- 
sinh scale ranging from to 90, with steps of 10 in M© pc"^ 
with CO contours at levels 0.8, 2, 4, 6 and 8 K kms/s from 
[Fukui et al. (2008)( Average intensity of the radiation field, {U), 
in rainbow arcsinh color scale ranging from to 20, with steps of 
4 in units of the local solar neighborhood intensity. Distribution 
of the PAH-to-dust mass fraction, fpAH, in units of the Galactic 
fraction of 0.046, as measured in the pixel based dust modeling 
with a rainbow color linear scale ranging from to 0.8 in steps of 
0.2 with the comparable fpAH from |Paradis et al. (2009)| in black 
contours levels of 0.007, 0.2 and 2. 
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